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We investigate the behaviour of stochastic differential equations, especially 
Burgers' eq., by means of Monte-Carlo-techniques. By analysis of the produced 
configurations, we show that direct and often intuitive insight into the funda- 
mentals of the solutions to the underlying equation, like shock wave formation, 
intermittency and chaotic dynamics, can be obtained. We also demonstrate 
that very natural constraints for the lattice parameters are sufficient to ensure 
stable calculations for unlimited numbers of Monte-Carlo-steps. 
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1. Introduction and Motivation 

Hydrodynamic turbulence remains a basically unsolved problem of modern 
physics. This is especially noticeable as the fundamentals seem to be fairly 
easy — the Navier-Stokes-equations for the velocity and pressure fields v 
and p, 

dtVa + VjsdpVa - l^'^^'^Va + ^daP ^ 0, (1) 

with the additional constraint 



daVa = 0, 



(2) 
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express the conservation of momentum in a classical, newtonian, incom- 
pressible fluid of viscosity v and density p. Laminar flows are reproduced 
very accurately; in the turbulent regime, it still is an open question how the 
universal characteristics of a flow, the scaling exponents of the structure 
functions Sp of order p, defined by 

Sp{x):={[Hr + x)-v{rW),^\x\^^. (3) 

can be extracted from the basic equations. Intermittency is reflected by 
exponents ^j, that differ from those expected by dimensional analysis. 
Monte-Carlo-simulations enable us to analyze turbulent flow patterns in 
detail, to gain direct insight into the formation of localized structures and 
their behaviour, and to measure observables like the scaling exponents. 
Instead of working with the full Navier-Stokes-eqs., we decided to elaborate 
the methods using Burgers' eq. 

dtVoi + VfjOpVa - vV'^Vct = 0, (4) 

which can be interpreted as the flow equation for a fully compressible fluid. 

A finite viscosity ly and energy dissipation e provide a dissipation length 
scale A corresponding to the Kolmogorov-scale in Navier-Stokes-turbulence: 




Besides being of interest of its own (e.g. in cosmology), working with eq. 
(4) has a number of technical advantages:* 

• Burgers' eq. is local, while the incompressibility condition acts as 
a nonlocal interaction in Navier-Stokes-turbulence. 

• The fundamental solutions to Burgers' eq. are well-known; in the 
limit of vanishing viscosity (Hopf-eq.), these form singular shocks as 
seen in fig. (1). The dissipation-scale provides an UV-regularization 
of the shock structures. 

• A huge variety of analytical methods have been applied to Burg- 
ers eq., giving results that can directly be compared to numerical 
measurements; and the origin of intermittency is well understood. 

2. Path Integral 

For the moment, we restrict our work to 1-dimensional Burgulence which 
already shows intermittent statistics. A random force / has to be introduced 

"For a more complete overview, see^*^ ^"^ . 
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Fig. 1. Typical solution of Burgers' eq. in the limit of vanishing viscosity. The graph 
shows v{x, t) as a function of x in periodic boundary conditions at constant t. 



into eq. (4) so that solutions are statistically homogenous in time: 

d d"^ 

B[v\ = dtv + v—v-v-^v = f, (6) 

where we model the stochastic force to be gaussian distributed with zero 
mean, energy injection rate e and correlation length A: 

(/(^,*))=0, (7) 

X-\x,x'-,t,t') = {f(x,t)f{x' ,t')) = e5{t-t')e^^(-''-^^-^\ . (8) 



A 

We then expect intermittent statistics to be found within the inertial sub- 
range A <C a; <C A. We can further identify a characteristic velocity at the 
injection scale, uq = (eA)^/^, and a Reynolds-number 

Re = (eAVz^^)'/^ (9) 

The path integral is introduced via the standard Martin-Siggia-Rose- 
formalism, giving for the generating functional*^ 



(10) 



Z[J] = J VvVfS{B[v]-f)exp (^-i j fxf + j Jv 

= J Vve-^^''''^, (11) 
with the action 

S[v;J]= j dxdtdx' dt' {B[v{x,t)\x{x,x' ■,t,t')'B[v{x' ,t')\ + J{x,t)v{x,t)) . 

(12) 



The functional determinant can be shown not to contribute for local theories, see 

_ _. HMPV 99 
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Beginning from an equivalent path integral, it has been shown that inter- 
mittent statistics of Burgers' eq. can be understood in terms of instanton 

solutions (I^FKL 97) 



3. Monte-Carlo-Simulations 

We discretizcd the above action onto a rectangular lattice with L sites in 
space-, and T sites in time-direction. Derivatives have been written in a 
symmetric (Stratanovich-) prescription. Lattice spacings will be denoted 
Aa; and At, respectively. We mainly used a heat bath algorithm on single 
nodes. 



3.1. Lattice parameters 

To indentify the lattice parameters with the constants of the continuum 
theory, we first notice that the viscosity has to be defined as"^ 

in which we define the arbitrary constant a = 1. The so-defined u gives us 
Re according to eq. (9). We further find that the dissipation length A is 
related to the correlation length A simply by 

In any practical application, it is sufficient to define the ry and Re, to chose 
L and A, and to calculate from that T and e. Stability considerations lead to 
further constraints for the lattice size, as will be explained in the following 
subsection. 



3.2. Stability 

As would be expected, the stability of the simulations over a large number 
of Monte-Carlo-steps is a big issue, due to the shock-like solutions of eq. (4). 
Indeed, if certain restrictions to the lattice parameters are not taken care 
of rigorously, the simulation terminates sooner or later due to divergencies. 
It is interesting to notice that the occurence of instabilities in our MC- 
simulations is related to the (non-trivial) existence of the dissipation scale 
A. We found that to obtain stable simulations, A has to be resolved on the 



''For example, this can be shown via the continuum Umit of the symmetric random walk, 
leading to the diffusion eq.. 



August 26, 2008 18:18 



WSPC - Proceedings Trim Size: 9in x 6in homeier 



5 

lattice: A > A.x. Unstable simultations occur otherwise — we observed that 
the overall energy of the configurations accumulate in the smallest scale Ax, 
causing the configuration to separate into two sub-configurations, of which 
one looks as the expected solution to the Hopf-eq., while the other grows 
beyond any limit, eventually breaking the simulation. 

As long as the dissipation length is resolved, the simulations are stable. Hav- 
ing performed several millions of MC-steps, no further instabilities occured. 
If length scales arc measured in units of the system size, this translates into 
a constraint involving the Reynolds-number: 

Re < AL, (15) 

which enforces huge lattices for high i?e as A < 1. 

3.3. Configurations 

We simulated systems of different sizes, from (L = 4) x (T = 16) to larger 
lattices of the same viscosity, as 8 x 64 and 16 x 256, and also different 
viscosities, like 16 x 16 or 64 x 32. The Reynolds-number ranged from Re = 
0.01 to Re = 100, not respecting the above conditions for stability. This 
may seem surprising, but we could see that interesting information could 
also be extracted from the physical sub-configurations. Stable runs have 
been obtained for Re = 1. 

Long runs of several millions of MC-steps, on a single node, are realistic for 
small lattices like 4 x 16 or 8 x 64 only. For larger lattices, a parallelized 
version of the code will be needed. 
We obtained the following results: 

• Thermalization and autocorrelation times are very long, up to the 

order of 10^^ MC-steps. 

• In the stable runs, after thermalization, the typical shock solutions 
of the Burgers' eq. form and can be observed moving and interact- 
ing with each other, see fig. (2). 

• In the unstable runs before occurence of the instability, configura- 
tions resemble the kink-solutions of the Hopf-eq.. 

• The disctinction between stable and unstable simulations can di- 
rectly be related to the existence of a dissipation length scale; whic;h 
is cither bigger (stable) or smaller (unstable simulations) than the 
lattice spacing. 
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Fig. 2. Calculated configuration v{x, t) of a (L = 32) X (T = 8)-lattice; the plot shows 
a slizc of constant t, periodic boundary conditions in x. Two shock-like structures are 
clearly visible- 



4. Summary and Outlook 

We have shown how to perform stable MC-simulations of stochastic partial 
differential eq., like Burgers' or the Hopf-eq.. The lattice versions of the 
theories can directly be identified with their continuum counterparts, and, 
as long as certain constraints on the lattice size are respected, unlimited 
numbers of configurations produced. Direct insight into the structures lead- 
ing to intermittency and, thus, multiscaling, can be obtained. Especially, 
wc want to point out that the existence of a dissipation length scale can be 
observed. 

As next steps, complete statistics will be made; especially the scaling expo- 
nents of the structure functions have to be compared with the analytic re- 
sults. Later, we will proceed by analyzing the incompressible Navier-Stokes- 
eqs.. 
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